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MODEE PREDICTIONS 
OF MAXIMUM THEORETICAL YIELDS OF AMINO AC IDS FOR 
GROWTH ON GLUCOSE AND ACETATE 



Maximum Theoretical Yield Maximum Theoretical Yield 

(mmol / per 10 mmol Glucose) (mmol^ per 10 mmoi Acetate) 





Palsson 
'93 


Modified 
Keasling '97 


Universa 
I Model 


% Increase! Pal ^ on 


Modified 
Keasling '97 


Universal 
Model 


/o 1 1 lOi trdb U 


Alanine 


20.00 


~ 20~. 00 


20.00 




3 93 


5,29 


" 5 29 


- 


Arginine 


7.74 


9.26 


10 07 


8.75 c 'o 


1 51 


2 43 


2 65 


9.05% 


Asparagine 


15.60 


18.18 


19 23 


5.77% 


3 24 


4 66 


4.91 


5.45% 


Aspartate 


18.20 


20 00 


20.00 




3 82 


5 29 


5 29 




Cysteine 


9.75 


1 1 49 


1 1 90 


3.57% 


1 81 


3 29 


3.42 


3.80° o 
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10.00 


13 33 


13 33 




2 68 


3 65 


3 65 




Glutamine 


10.00 


13 33 


13 33 




2 50 


3 46 


3 46 




Glycine 
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35 33 


35 33 




3 94 


9 00 


9 00 




Histidine 


7.30 


9.77 
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0 23 c o 


1 37 
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7.34 
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0 91% 


1 44 
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5.87 




Threonine 


12 30 


15 00 


15.00 




2 50 


3 91 
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10 00 
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1 96 
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Palsson '93: E. coli model proposed by Palsson ( 1 993) 

Modified Keasling '97: Modified Keasling ( 1997) E. coli model as described in text 

Universal Model: Modified Keasling (1997) E. coli model augmented with non-/:. coli reactions 

compiled by the Kyoto Encyclopedia of Genes and Genomes 
% Increase: Betw een the modified Keasling ( 1997) model and the Universal model 



Fig. 9 
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MODIFICATIONS TO THE PRAM ANIK AND KHASLING MODEL* 



Enzymes 

Reactions assumed irreversible 

Phosphofructokinase 

Citrate Synthase 

2-Ketoglutarate Dehydrogenase 

PRSCAIM Synthetase 

Glycerol Kinase 
Reactions removed from model 

Unknown Pathway 

Cystathionase 

Sulfotransferase 
Reactions modified 

Fructose-1 ,6-bisphosphate Aldolase 

Isocitrate Dehydrogenase 

Succinate Thiokinase 

Prephenate Dehydrogenase 

Hoi Dehydrogenase 

RCAIM Synthetase 

GTP Cyclohydrolase 

3,4-Dihydroxy-2-Butanone-4-Phosphate 
Synthase 

H2Neopterin Triphosphate 
Pyrophosphatase 

CoA Synthase 



Reactions 



Pi 



Fructose-1 ,6-bisphosphate -> Fruc tose-6-phosphate 
Acetyl-CoA + Oxaloacetate -> CoA + Citrate 
2-Ketoglutarate + NAD + CoA ~> Succinyl-CoA + C02 + NADH 
RCAIM + ATP + Aspartate ADP + Pi + PRSCAIM 
Glycerol + ATP -> Glycerol-3-phosphate + ADP 

S'-methylthioadenosine -> Adenosine + Methionine 
Homocysteine + Adenosine <- s-Adenosyl-homocystine 
Adenosine-3,5-diphosphate + sulfite <---> 3-Phosphoadenylylsulfate 

Fructose-1 ,6-bisphosphate -> Fructose-6-phosphate + Pi 
Isocitrate + NADP < — > C02 + NADPH + 2-Ketoglutarate 
Succinyl-CoA + ADP + Pi < > ATP + CoA + Succinate 
Prephenate + NAD -> C02 + NADH + para-Hydroxy phenyl pyruvate 
Histidinol + 3 NAD -> 3 NADH + Histidine 

AIR + C02 + ATP -> 5-p-Ribosyl-4-carboxy-5-aminoimidazole + ADP 
+ Pi 

GTP ~> D6RP5P + Formate + Ppi 
Ribulose-5-phosphate — > DB4P + Formate 

AHTD - > PPi + Pi + DHP 

OIVAL + METTHF + NADPH + ALA + CTP + 4 ATP + CYS -> 
THF + NADP + AMP + 2 PPi + 2 ADP + CQ2 + CoA + CDP 



MODIFICATIONS BASED ON INFORMATION BY KARP (1999) 



Fig. 12 
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''' ' Same yetie responsible for two intracellular reactions 
I,K Same yene responsible tor [wo intracellular reactions 
' No gene has been assigned to these intracellular reactions 



Fig. 13 
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GENES SELECTED FOR REMOVAL BY KNOCKOUT STUDY 



Enzymes 



Genes 



Reactions 





AP e 


/-\L.t? La it; r\ 1 1 I dot! 




OLJr rXlllabfcJ 


1 l(J r\ 


CMP Kinase 


ndk 


F0F1 -ATPase 


unc 


Formate THF Ligase 


FTL C 


Fumarase 


fumAB 


Glyceraldehyde Kinase 


GK d 


Glycine Cleavage System 


gcvHTP 


Malate Dehydrogenase 


mdh 


Methenyl THF Cyclohydrolase 


folD f 


Methylene THF Dehydrogenase 


folD 9 


NADH Dehydrogenase I 


ndh 


PEP Synthase 


pps 


Phosphatidate Phosphatase 


dgkA 


Phosphotransacetylase 


pta 


Pyrophosphatase 


ppa 


Succinate Dehydrogenase 


sdhABCD 


Succinate Thiokinase 


sucCD 



35ADP -> AMP + Pi 
AC + ATP -> ACTP + ADP 
CDP + ATP -■> CTP + ADP 
CMP + ATP -> CDP + ADP 
ADP + Pi + H ext -> ATP 

THF + FORMATE + ATP ADP + Pi + FTHF 
FUM -> MAL 

GLAL + ATP ADP + T3P1 

GLY + THF + NAD > METTHF + NADH + C02 + NH3 
MAL + NAD -> NADH + OA 
METHF -> FTHF 

METTHF + NADP -> METHF + NADPH 

NADH + Q > NAD + QH2 + 4 H ext 

PYR + ATP -> PEP + AMP + Pi 

DGR + Pi •-» PA 

ACTP + COA -> ACCOA + Pi 

PPi -> 2 Pi 

SUCC + FAD -> FADH2 + FUM 

SUCCOA + GDP + Pi -> GTP + COA + SUCC 



a,b Same gene responsible for two intracellular reactions 
f,g Same gene responsible for two intracellular reactions 

c,d,e No gene has been assigned to these intracellular reactions 



Fig. 14 
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MODEL SELECTIONS OF ENZYMATIC REACTIONS THAT WILL 
ENHANCE I I IE AMINO AC ID PRODUCTION CAPABILITIES OF 
ESCHERICHIA COLI 



Amino Acid Substrate 
Arginine Glucose: 

Acetate: 

A Glucose/ 
Asparag.ne Acetate; 

~ . . Glucose/ 
C * Steme Ace.ate: 
Histidine Glucose: 

Acetate: 



Isoleucine Glucose: 
Methionine Glucose: 
Acetate: 



Tryptophan Glucose: 
Acetate: 



EC# 

2.7.1.90 
2.7.2.2 
2.7.2.2 
2.7.2.12 

•...3.1.4 

2.7.7.5 

1.4.1.10 
2.7.1.90 
1.4.1.10 

4.1.1.38 

many 

2.7.7.5 

1.4.1.10 

2.7.7.5 

2.7.9.1 

4.1.1.38 

2.7.1.90 

2.7.9.1 

2.7.9.1 

4.1.1.38 



Enyzme 



6-Phosphofructokinase (pyrophosphate) 

Carbamate kinase 

Carbamate kinase 

Acetate kinase (pyrophosphate) 

Aspartate— ammonia ligase (ADP- 
fornung) 

Sulfate adenylyltransferase (ADP) 

Glycine dehydrogenase 
6-Phosphofructokinase (pyrophosphate) 
Glycine dehydrogenase 
Phosphoenolpyruvate carboxykinase 



(pyrophosphate) 

Sulfate adenylyltransferase (ADP) 
Glycine dehydrogenase 
Sulfate adenylyltransferase (ADP) 
Pyruvate, phosphate dikinase 

Phosphoenolpyruvate carboxykinase 
(pyrophosphate) 

6-Phosphofructokinase (pyrophosphate) 
Pyruvate, phosphate dikinase 
Pyruvate, phosphate dikinase 

Phosphoenolpyruvate carboxykinase 
(pyrophosphate) 



Reaction Catalyzed 

Fructose-6-P + PPi > Fru< :tose-1 ,6-Bisphosphate + Pi 
ATP + NH3 + C02 > ADP + Carbamoyl Phosphate 
ATP + NH 3 + C02 > ADP + Carbamoyl Phosphate 
Acetate + PPi > Pi + Acet/I-Phosphate 

ATP + NH3 + L-Aspartate > Pi + ADP + L-Asparagme 

Sulfate + ADP > Pi + Adenylyl-Sulfate 

NAD + glycine > glyoxylate + NADH + NH3 
Fructose-6-P + PPi > Fructose-1 ,6-Bisphosphate + Pi 
NAD + glycine > glyoylate + NADH + NH3 

PPi + Oxaloacetate -> CQ2 + Pi + PEP 



Sulfate + ADP - > Pi + Adenylyl-Sulfate 
NAD + glycine > glyo>ylate + NADH + NH3 
Sulfate + ADP -> Pi + Adenylyl-Sulfate 
Pyruvate + Pi + ATP -> AMP + PPi + PEP 

PPi + Oxaloacetate > C02 + Pi + PEP 

Fructose-6-P + Ppi -> Fructose-1 ,6-Bisphosphate + Pi 
Pyruvate + Pi + ATP --> AMP + PPi + PEP 
Pyruvate + Pi + ATP -> AMP + PPi + PEP 

PPi + Oxaloacetate - > C02 + Pi + PEP 



Fig. 15 
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Target % of Maximum Growth Rate 



Fig. 21 
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EVOLUTION OF MINIMAL REACTION SETS FOR C ASE (I) 
UNDER DECREASING GROWTH REQUIREMENTS. 



Target % 
.Maximum 
(; row th Rate 


Minimal 
Reaction Set 
(# Reactions) 


Key Features 


100% 
90% 


234 


The glycolysis, tricarboxylic acid cycle, and pentose phosphate pathways are all 
operating in their forward directions, optimally generating the energy cofactors ATP, 
NADH, and NADPH required for cell growth. All available glucose is oxidized into 
the cell's only secreted byproduct, carbon dioxide. 


229 


The fluxes through two TCA cycle reactions 2-ketogiutarate dehydrogenase and 
succinate dehydrogenase are zero while succmyl-CoA synthetase operates in its 
reverse direction suggesting a less demanding energetic state under the sub- 
maximal growth demands. Acetate is now secreted as a byproduct along with 
carbon dioxide. 


80% 


228 


Fluxes through two additional TCA cycle reactions, furnarase and malate 
dehydrogenase, are eliminated while a reaction secreting fumarate is added. 


70% 


226 


The pentose phosphate pathway operates solely for nucleotide biosynthesis with the 
reaction fluxes through ribulose phosphate 3-epimerase, transketolase I, 
transketolase II, and transaldolase B all operating m reverse. Fluxes through 
glucose-6-phosphate dehydrogenase, lactonase, and 6-phosphogluconate 
dehydrogenase are absent in this case, replaced by pyridine nucleotide 
transhydrogenase which meets the cellular NADPH needs. In addition, formate is 
now secreted along with acetate, fumarate, and carbon dioxide. 


60%, 50%, 
40% 


225 


Acetate is no longer secreted as a metabolic byproduct, but is converted to acetyl- 
CoA by acetyl-CoA synthetase. 


30%, 20%, 
10%, 1% 


224 


Three glycolytic reactions, phosphoglycerate mutase, enolase, and pyruvate kinase 
are eliminated, but both serine deaminase and phosphoenolpyruvate synthase are 
added to supply the cell with phosphoenolpyruvate. 



Fig. 22 
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METABOLITES UPTAKEN OR SECRE TED 
RATE ON AN OPTIMALLY ENGINEERED 
U DENOTES METABOLITE 
S DENOTES METABOLITE 
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Perce n tage of 100% Biomass Generation Required 



Metabolite 


100% 


99.5% 


99% 


98% 


! 97% 


] " ■ 
! 96% i 95% 


! 90% ! 85% ! 80% '70% '60% 


10% 


Acetate 

Acetaldehyde 
Adenine 







1 

; u 


; . I j' " j" " " ; 

u ! u ! u ! u i u 


s 

: U ! U ! 


s 

! U 

\ 

! u _ 

I 

u 


Adenosine 




_j _ 












Alanine 











T 




u 


u 

u i u 


I u : 


Arginine 


u 


u u 


r 

! U 


u 


u ; u 


; u 


! U 


Asparagine 


















u ! u _) 


u 


__u 
u 
s 
u 


Aspartate 




s 














U ! U i U 


_ u _ 


Carbon dioxide 


s 


s 


s 


s 


s 


s 


s 


sjs 


; s 


: s 


Cysteine 


u 


u 


u 


u 


u 


u 


u u ! u ! u 


U ! U 


D-Alanine 








I— 




U i u I 


r 

1 u 


_u__ 


Thymidine 




u 


u 


u 


__u 

u 


u 


u r u 


1 u 1 u 
* u 


r - 1 " 
U ! U 


_u 


Ethanol 


u 


u 


u 


u 


u 


r u < u 




u 




Glycerol 














u 








u 






Glycerol-3-phosphate 


u 


u 


u 


u 


u 


u 


U 


u 


u 




u 


u 


Glutamine 




















u 


_u_ 


_u_ 


_u_ 


Glutamate 






















s 


u 


u 


Glycine 






1 

1 




u 


u 


u ! u ; u 


u 


~ — - ■ 

u 


, . 

u 


Guanine 








u 


u 


u 


U ! 


u 


r~ — 

u 








Guanosine 
















U 






u 


u 


u 


Histidine 




u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Isoleucine 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 
u 


_u _ 
u 


Leucine 














u 


u 


u 


u 


u 


Lysine 


u 


u 


u 


u 


u 


u 


u 


U 


u 


u 


u 


u 


u 


Meso-diaminopimelate 




u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Methionine 


u 


u 


u 


u 


u 


u 


u 


u 


U 


u 


u 


u 


u 


Mannitol 
























u 


u 


Ammonia 


u 


u 


u 


u 


u 


u 


u 


u 












Oxygen 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Phenylalanine 






u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Phosphate 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Proline 










u 


u 


u 


u 


u 


u 


u 


u 


u 


Putrescine 


u 


u 


u 


u 


u 


u 


u 


u 


U 


u 


u 


u 


u 


Pyruvate 




















u 


u 


u 


u 


Ribose 
























u 


u 


Serine 
















u 


U 


u 


u 


u 


u 


Spermidine 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Threonine 




u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


u 


Tryptophan 




u 


u 


u 


u 


u 


u 


U | u 


u 


u 


u 


u 


Tyrosine 






u 


u 


u 


u 


u 


U U 


u 


u 


u 


u 


Uracil 












u 


u 


U u 






u 




Uridine 




















U | 




Valine 














u i U | u 


u 


Si 




u 


# Metabolites Uptaken 


12 


17 


19 


21 


22 


24 


26 


28 


29 31 


29 


34 


34 



Fig. 23 




100% 99% 98% 97% 96% 95% 94% 93% 92% 91% 90% 



Target % of Maximum Growth Rate 



Fig. 24 



^^ttC 

m 



Applicant: MARANAS, Costas D., et al. 
attorney Docket No. P05468US1 

id: January 10, 2002 
w le: METHOD & SYSTEM FOR MODELING 
CELLULAR METABOLISM 
Sheet 23 of 26 




Fig. 25 



^1 



Applicant: MARANAS, Costas D., et al. 
ttorney Docket No. P05468US1 

led: January 10, 2002 
tie: METHOD & SYSTEM FOR MODELING 
CELLULAR METABOLISM 
Sheet 24 of 26 



EVOLUTION OF MINIMAL REACTION SITS FOR C ASE (II) 
UNDER DECREASING GROWTH REQUIREMENTS. 



Target % 
Maximum 
Growth Kate 


Minimal 
Reaction Set 
(# Reactions) 


Key Features 


100% 


201 


The organic material transported into the cell includes ethanol and glycerol-3- 
phosphate which fuel glycolysis, the TCA cycle, and PPP The flux directions of the 
glycolysis pathway are split with all reaction fluxes preceding glyceraldehyde-3- 
phosphate (G3P) dehydrogenase operating in reverse, and all flu<es following and 
including G3P dehydrogenase operate in their forward directions. Putrescine, 
spermidine, and five amino acids are transported into the network eliminating the 
need for biosynthetic pathways for these components. 


90% 


132 


While the PPP and TCA cycle reactions are still functional, the network no longer 
utilizes the five glycolytic reactions from glyceraldehyde-3-phosphate 
dehydrogenase to pyruvate kinase. Consequently, the TCA cycle is completely 
fueled by imported ethanol and acetate rather than flux from the glycolysis pathway. 


80% 


125 


This network tolerates the complete elimination of the TCA cycle and glyoxylate 
shunt. As a result, the function of the pentose phosphate pathway reactions is no 
longer restricted to nucleotide biosynthesis, but now includes the formation of 
cellular NADPH. Most of this NADPH is subsequently converted to NADH by 
pyridine nucleotide transhydrogenase to replace the cellular reducing power lost 
from the inactivity of the TCA cycle. 


70% 


124 


A slightly less efficient set of internal metabolic reactions enables the growth 
demands to be met with the importation of one less metabolite (i.e. one less 
transport reaction) than its 80% counterpart. 


60% 
50%, 40% 
30%, 20% 


123 


Neither the TCA cycle nor PPP are utilized for reducing power. Most of the cellular 
reducing capabilities are now generated from the uptake of ethanol and its 
subsequent conversion into acetyl-CoA. 


10%, 1% 


122 1 


This minimal network is comprised mostly of cell envelope and membrane lipid 
biosynthetic reactions, along with a number of transport and salvage pathway 
reactions. Here, the three core metabolic routes, glycolysis, the TCA cycle, and the 
pentose phosphate pathway are almost completely dismantled with only one 
glycolytic and 4 PPP reactions remaining. 
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FUNCTIONAL CLASSIFICATION OF MINIMAL 

NETWORK REACTIONS FOR GROWTH ON AN OPTIMALLY 

ENGINEERED MEDIUM. 



Functional Classification # rxns 


ALA Isomerization 1 
Alternative Carbon Source 7 


Anaplerotic Reactions 


1 


Cell Envelope 
Biosynthesis 


29 


EMP Pathway 
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Membrane Lipid 
Biosynthesis 
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Pentose Phosphate 
Pathway 
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COMPARISON OF MINIMAL METABOLIC GENE/REACTION 
SETS BASED ON FUNCTIONAL CLASSIFICATION 



Metabolic Function 


Essential 
Gene Set 
Ref. (2) 


Minimal 

Gene Set 
Ref. (5) 


Minimal 
Reaction 
Set 


# Genes 


# Genes # Reactions 


Amino acid biosynthesis 0 

Biosynthesis of cofactors, prosthetic 

groups, and carriers i 4 


<> j 1 

i 

3 ! 0 


Cell envelope 2 


11 1 29 


Central intermediary metabolism 


7 


7 


1 


Energy metabolism 


31 


32 


r 21 


Fatty acid and phospholipid 
metabolism 


5 


7 


16 


Purines, pyrimidines, nucleosides, 
and nucleotides 


17 


14 


18 


Transport and binding proteins 


17 
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